NATIONAL ADVISORY COMMITTEE 

FOR AERONAUTICS 

/ 

TECHNICAL NOTE 2499 

LAMINAR FRICTION AND HEAT TRANSFER AT MACH 

NUMBERS FROM 1 TO 10 

By E. B. Klunker and F. Edward McLean 

Langley Aeronautical Laboratory 
Langley Field, Va. 





LIBRARY KAFB, NM 


TECH UBRAHY KAFB, MM 


1 


llllllll 

I 

I 

ill 

III 

ill 


]□ 


D1 



HATIONAL ADVISORY COM4ITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2^99 

LAMINAR FRICTION AND HEAT TRANSFER AT MACH 
NUMBERS FROM 1 TO 10 
By E. B. Klunker and F. Edward McLean 

SUMMARY 


A method of solving the laminar boxmdary- layer equations for com- 
pressible flowj in the absence of a pressure gradient, is. developed 
without imposing restrictions on the thermal properties of the fluid 
niedium. Velocity and temperature profiles and boundary- layer charac- 
teristics have been con^juted for Mach numbers from 1 to 10, utilizing 
experimental values of the heat capacity, viscosity, and conductivity. 
The analysis shows that an effective temperature, which is a function 
of the surface temperature and stream conditions and is similar to the 
recovery temperature, arises naturally and is the proper reference tem- 
perature to be Tised in heat-transfer calculations. The effective tem- 
perataire and the recovery tenqjerature become identical for the condition 
of zero heat transfer. The recovery factor and the analogous effective- 
temperature ftinctlon decrease substantially with increasing values of 
Mach number. Thus, the high-speed aerodynamic heating problem is not 
as severe as is indicated by simpler theories. 


INTRODUCTION 


The large extent of laminar flow expected on missiles and aircreift 
operating at high speeds and altitudes makes an accurate knowledge of 
the laminar-flow frictional effects and heat-transfer characteristics 
desirable for practical design considerations. The calculation of the 
bOTondary- layer characteristics at high speeds is complicated by the 
necessity of considering simultaneously the compressibility of the fluid 
medixmi and the variation of thermal properties with temperature. As a 
consequence of the conipresslble nature of the flow, the tenperatiire 
variation throughout the boundary layer is large and the variation of 
the thermal properties of the fluid is significant. At a stream Mach 
number of 5> for eianqjle, the ratio of maximum boundary- layer teu 5 )eratTire 
(absolute) to stream ten 5 )erature is of the order of 5 to 1 for the 
condition of zero heat transfer at the surface. Moreover, for these 
conditions, experiment shows that the viscosity and conductivity 
increase by factors of approximately 3 and k to 1, respectively, over 
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their values at the edge of the houndary layer, while the heat capacity 
increases hy approximately l6 percent and the Prandtl number decreases 
more than 10 percent. Thus, since the tengierature variation within 
the houndary layer at hi^ supersonic speeds is large, assumptions or 
restrictions regarding the thermal properties can conceivably have a 
significant influence on the calculated characteristics of the houndary 
layer. 

The calculation of the characteristics of the laminar hoiondary 
layer at high speeds has been considered by several investigators 
(references 1 to 4, for exan^jle). In these and similar analyses, how- 
ever, certain restrictions are Imposed on the thermal properties of the 
fluid. Two assun5)tions are generally made: First, the heat capacity 

and the Prandtl number are considered to be independent of temperatiore; 
and, second, some simple law for the variation of viscosity (or con- 
ductivity) with temperature is specified. The constancy of the heat 
capacity and Prandtl number implies that the viscosity and conductivity 
vary with temperature in the same manner. For air, these conditions 
are nearly satisfied for small temperature variations. The viscosity- 

temperature relation is frequently taken as the equation — = 

^"l 

where ou is a constant, a linear relation between viscosity and tem- 
perature, or the Sutherland equation. Of these three, the Sutherland 
equation especially gives good results for the viscosity over a wide 
temperature range, although it is not as satisfactory for the conduc- 
tivity. The choice of a linear viscosity- temperature relation (refer- 
ence 4, for example) or a Prandtl ntmiber of unity (reference 1, for 
example) leads to an essential slmplifipation, since in each case only 
one differential equation need be solved. Although these vario\zs 
restrictions on the thermal characteristics appear justified for small 
temperature variations, they are not justified a priori for the calcu- 
lation of boundary- layer characteristics at high Mach numbers where the 
temperature variations are large. 

The present investigation is concerned with the calculation of the 
characteristics of the laminar boimdary layer in the absence of a pres- 
sure gradient at supersonic speeds. The purpose of the investigation 
is twofold: first, to determine the characteristics of laminar boundary 

layers at supersonic Mach numbers and, second, to provide information 
which may be \ised to evaluate restrictions frequently imposed on the 
thermal characteristics. The method developed herein places no 
restrictions on the thermal properties of the fluid. Experimental 
values for the variation of heat capacity, conductivity, viscosity, and 
Prandtl number have been employed in the calculations in order to obtain 
the flow properties as accurately as possible. The method can be adapted 
to the calculation of boimdary layers with suction or blowing. 
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Computations have been made for a surface ten 5 )erature equal to 
the stream temperature at Mach numbers from 1 to 10 and for an Insulated 
surface for Mach numbers from 1 to 5. In addition, calculations have 
been made at a Mach n\jmber of 7 for ratios of surface teraperatiare to 
stream temperature of 2 and 4. The stream temperature was taken 
as -67° P, a value that coiresponds to the isothermal region of the 
MCA standard atmosphere. 

The method of solution is essentially one of successive approxi- 
mation. Ten iterations were required in most cases to obtain a dif- 
ference of less than 1 percent, for both the shear stress and the 
heat-transfer rate at the surface, between successive approximations. 
Even with the use of automatic computing equipment (the Bell Telephone 
Laboratories X- 667 ^ relay conqruter at the Langley Laboratory was used) 
appreciable time is required for the solution. As a consequence, the 
number of computations con^jleted thiis far is inadequate for a critical 
comparison of calculated boiindary- layer characteristics with other 
analyses . 


SYMBOIS 


x,y coordinates parallel and normal to stream direction 

ijf stream function 


^PlUlUqX 



u,v components of velocity along x- and y-axes 

Uq stream velocity 

M Mach number 

T temperature 

d temperature function 

p density 
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V- 

V 

k 

q 

li 

R 

Pr 

Nu 

T 

D 

Cf 




keat capacity at constant pressure 

ratio of heat capacities at constant presstire and 
constant volume 

coefficient of viscosity 

kinematic viscosity 

thermal conductivity 

local heat-transfer rate 

local heat-transfer coefficient 

Reynolds number 

Prandtl number 

Russelt n-umher 

shear stress 

skin friction 

skin- friction coefficient 

integrating factors 


f = p*\i*n* 

C-p* 

Pr 

function defined hy equation (l3) 
F 2 function defined hy equation (16) 

g function defined hy equation (16) 

cp function defined hy equation (l4) 

fimction defined hy equation (l7) 


<I> 
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Subscripts: 

1 free stream or at edge of boundary layer 

e effective 

r recovery 

s surface 

Superscript: 

* dimensionless quantity based on stream conditions 

AMLYSIS 


Basic equations .- With the x-axis taken in the free-stream direc- 
tion^ the two-dimensional compressible-flow boundary- layer equations 
for steady motion in the absence of a pressure gradient may be written 
as follows; 

Momentum equation: 

pv 2a = kfp 

3x Sy . ^y/ 

Continuity equation: 


(pu) + (pv) = 0 
ox oy 


Energy eqxaation: , 


( 1 ) 


( 2 ) 


( 3 ) 
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where u and v are the velocity components in the direction of the 
X- and y-axes^ respectively, p is the density, T is the absolute 
temperature, Cp is the heat capacity at constant pressure, p, is the 

coefficient of viscosity, k is the thermal conductivity, and the 
thermal properties are functions of the temperature. 

The set of partial-differential equations (l) to (3) may be 
reduced to two ordinary differential equations for the particular cases 
of constant tenrperature or zero heat transfer at the surface y = 0. 

In order to accomplish this reduction it is convenient first to intro- 
duce the stream function as a new independent variable. The velocity 
components are related to the stream function by the equations: 




Then by taking x and ijr as the independent variables, equations (l) 
and (3) may be transformed into 



and eqiiatlon (2) is satisfied Identically by virtue of equation (4). 

It is found that surfaces of constant fjp- are both isotherms 
and surfaces of constant velocity. Thus, with the choice of an inde- 
pendent variable proportional to the equations of motion and 

energy niay be reduced to ordinary differential equations for boundary 
conditions of constant temperature or zero heat transfer at the sur- 

i}f 

face y = 0. By introduction of the independent variable ^ = - 

^PlPlUiX 

and the nondimens ional temperature f line t ion ■fl, equations (5) and (6) 
become 


■ 5y 


" 'Sx 
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da dW 2 


( 7 ) 


. p^« II . -2pw(ff)' (8) 


where u* = 


u • .p, p T 


'k* — 


, T* = — , \i* = cip* = -^, Pr = 

’ fT_ ’ M - ^ Ct, ’ V ' 


,U. 


Pi 


T. 


Pi 


Pi 


•a = 


T* - 1 T* - 1 


31- Ij3m,2 

2cp^Ti 2 1 


The symbol y represents the ratio of heat capacities at constant 
pressure and constant voltnne, is the stream velocity, and the 

subscript 1 refers to values at the edge of the boundary layer or 
values in the free stream. 

Von Mises (reference 5) employed ilf as an independent variable 
in the boundary- layer equations for incompressible flow and many authors 
since have utilized this transformation for various compressible-flow 
problems. The transformations employed herein follow those of 
Von Karm& and Tsien (reference l). 

The solution of the boimdary- layer equations for conipressible flow 
may be simplified in tTO cases where the heat capacity and Prandtl 
number are takeji as constants. For a Prandtl nTomber of unity, the 
temperature and velocity have a simple parabolic relation and only the 
differential equation of motion need be solved to determine the flow 
throughout the boundary layer (reference 6). For a linear vlscosity- 
tenrperature relation the equation of motion is Independent of the 
energy equations and may be reduced to the form occurring in incompres- 
sible flow, (See reference k, for example.) Thus only the energy 
equation need be solved, since the incompressible-flow solutions are 
knovm. Where no restrictions are placed on "the heat capacity, viscosity, 
conductivity, and Prandtl number, both eq\iations (7) and (8) must be 
solved to determine the flow properties in the boundary layer. 

Although the nonlinear equations (7) and (8) do not yield to a 
direct analytical solution, they can be solved numerically by a method 
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of successive approximations. Since Cp*, [o.*, p*, and Pr are 

functions only of the temperature, and thus functions of the 
coefficients of the derivatives in equations (7) and (8) may he con- 
sidered as known functions if some initial approximate solution is 
given. Then, ulth the quantities 


f(^) = p*p.*u* 


and 


Ut) = ^ p*p*u* 

considered as known functions of and with the integrating factors 



the equation of motion (7) and the energy equation (8) may he written as 



Solution of houndary value problem .- Two problems are considered 

herein. The first is to determine the flow throughout the laminar 
boundary layer subject to the condition of a constant surface tem- 
perature, and the second is to determine the flow for the condition 
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of zero heat transfer at the siirface. The boundary conditions to be 
satisfied for these two problems are: 

For constant surface temperature: 

At the STorface {^ = O) , u* = 0 and i3 = ds = Constant 

At infinity ( ^ — > <») , u* = 1 and d = 0 

For zero heat transfer at the surface: 

At the svLTface (^ = 0), u* = 0 and 

= Constant ( p*u* — ] = 0 

V 

At infinity u* = 1 and d = 0 

The subscript s denotes the values at the surface. 

Since the boundary conditions on u* are the same for both problems, 
the solutions are of the same form. 

The first integral of equation (9) is 



of = 


and a second integration gives the velocity as 


u* = Fi ^ (12) 

Jo of 

From the boundary condition at infinity the function F^ is found to be 



The velocity distribution is thus determined from equations (ll) to (l3) 
for a given surface temperature (or for zero heat transfer at the 
surface) and for given stream conditions. 
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The temperature distrihution corresponding to a constant value of 
the STirface temperature is obtained hy integrating equation (lO). The 
first integral is 


= ( 1 ^) 

■where 

After a second integration the tenq)erature function is found to he 



From the hotindary condition at infinity, the function F 2 is deter- 
mined as 

Fg = g( 3 e - -eg) (16) 


•where 


s(»l.,Tl,Ts) 



The temperature distribution corresponding to zero heat treinsfer 
at the surface ^ = 0 is found in a similar manner. The first integral 
of equation (lO) satisfying the conditions of zero heat transfer is 
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where 


and after a second integration the temperature function is foimd to he 
For t = 0, ^ becomes the recovery factor 



Equations (12) and (l5) give the velocity and temperature distri- 
butions as functions of ^ corresponding to a specified surface tem- 
perature^ whereas equations (l2) and (l8) give the velocity and tem- 
peratvire distributions as functions of ^ corresponding to zero heat 
transfer at the solid surface. These solutions may be expressed in 


terns of the more familiar independent variable 



(where 


V 


is the kinematic viscosity) by means of the relation 


^ = p*u* 

dT] 



found from the expressions defining q and ^ and from equation (4), 

Comments on method of solution .- The solution of the equation of 
motion and the energy equation described herein imposes no restrictions 
on the variation of the material properties of the fltiid - heat capacity, 
viscosity, and conductivity - with temperature. Since the method of 
solution is essentially one of numerical Integration, experimentally 
determined values of these peirameters or values found from analytical 
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expressions may "be used witli equal facility. The eq\iations developed 
herein may he slnqilified somewhat, with an attendant reduction in 
conrputational labor, for the particuleir case in which the heat capacity 
the Prandtl number Pr are taken as constants. With these restric- 
tions the function X reduces to f/Pr and the integrating factor P 

Pt* 

reduces to a . 

The solution of the equations for the velocity and temperature dis- 
tributions throughout the boundary layer is a step-by-step method of 
successive approximations. Thus, initial approximate solutions for the 
velocity and temperature distributions are required. These initial 
solutions may conveniently be taken as those for Pr = 1 and p*|i* = 1 
or, better still, the solutions given in reference 4. These values are 
\ased to compute the quantities f, X, a, and p as functions of 
and the’ velocity and tenperature distributions are found by numerical 
integration. With these new values for the velocity and temperature 
distribution, the functions f, X, a, and p are recalculated and 
the second approximation to the velocity and temperature distributions 
is determined. The process is continued until the desired accuracy is 
obtained. This method of solution is essentially the same as the one 
enployed for evaluating the velocity distribution in reference 1. 


RESUITS ATTD DISCUSSION 


The method described in the preceding section has been employed to 
calculate a number of velocity and temperature profiles and bOTondary- 
layer characteristics at supersonic Mach numbers. Faired experimental 
values of cp, p, and k primarily from reference 7 were used in the 
conputatlons in order to obtain reliable results at the high Mach numbers 
considered. The free-stream tenperature Tq was taken as 392.7 Rankine 

(-67° F), the value at the Isothermal level of the standard atmosphere, 
and the corresponding values of the thermal properties were taken as; 

Cp^ = 7.718 Btu/( slug) (deg) 

Pq = 3.058 X 10-T slug/(ft) (sec) 

kp = 3.227 X 10"^ Btu/(ft)(sec)(deg) 

Prq = 0.73 


The caJ.c\alated velocity and temperature profiles for Tg* = 1 at 
Mach numbers from 1 to 10 are presented in figure 1, and some results at 
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a constant Mach number of 7 for values of Ts* =- 1^ 2, and 4 are given 
in figure 2. Profiles for the condition of zero heat transfer at the 
wall are shown in figure 3. Althoijgh the calculated temperatures are 
very high at high Mach numbers, it is significant that these values 
are considerably lower than those predicted by simpler theories. Also 
of interest is the fact that the velocity profiles for the condition 
of zero heat transfer are nearly linear, whereas the profiles for a 
constant surface temperatirre have appreciable ciurvature. The velocity 
and temperature approach their stream values at smaller values of t] 
for high Mach numbers than those predicted by siu^oler calcxilations; 
thus, a thinner boundary layer is indicated. 

The frictional and heat-transfer characteristics may be evaluated 
with the aid of the equations developed in the analysis. The shear 
stress at the wall is given by 


Ts 



Pi V / 

V 


du*\ 

A=0 


where R is the Reynolds nimiber From the equation defining a 

and from equation (ll), 



and the local skin friction may be written as 


Ts = 



( 19 ) 


The friction drag may be found by integration as 


D =/ 

J G 


X Pl^i^ 

Tg dx = 2 — — xFi 


w 
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and the average skln-frletlon coefficient then is given by 


D ^ 


( 20 ) 


The skin- friction coefficient found from equation (20) is shovn in 
figure 4 for the conditions of constant surface tempe'ratture and zero 
heat transfer. This figure shows that the skin friction is significantly 
influenced by both the Mach number and the thermal state at the s^lrface. 
The skin friction decreases with increasing Mach number, and at a given 
Mach number the skin friction increases with decreasing values of Tg*. 

At a Mach number of 5 for example, the skin friction decreases from the 
Mach number 1 value by approximately l8 percent for zero heat transfer, 
and for a value of Tg* = 1 the decrease is approximately 7 percent. 


The convective heat-transfer rate at the surface per uni t area (q) 
is given by the expression 






and from equations 


(l4) and (l6), since 


p = 1 and cp = 0 for 5 “ 0, 


(■ 






s 




(^e - 


2c 


Pi 


The local heat-transfer rate then may be written as 


q = -gc 


Pi 





( 21 ) 


The heat-transfer rate is often expressed as the product of a 
heat-transfer coefficient h and a temperature increment. The tem- 
perature increment Tg - Tg arises natturally in equation (21) eind 
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with this temperature difference an effective heat-transfer coefficient 
may he defined as 




PimUi 


( 22 ) 


The usual ten^ierature increment used for defining the heat- transfer 
coefficient is the difference between the recovery temperature and the 
surface temperature. With this temperature increment the heat-transfer 
coefficient becomes 


h = 


-gCp. 


1 


^e* - Ts* 


Te* - Tg* 

- Ts* 


(23) 


The relative magnitudes of the temperatiores T^, and Tg may he deter- 
mined from the recovei^ factor -3^, and the effective-temperature 
function 3g given in figure 5- Two points are significant: First, 

the recovery factor is considerably smaller than the effective-temperature 
fxmction; second, both ■flj. and 3g decrease with increasing Mach number. 

The value of the recovery factor is usually considered to be approxi- 
mately 0.85 for laminar flow. Figure 5 shows that for a Mach number 
of 5 the recovery factor is 0.766 and the effective-temperature function 
for Tg* =1 is 0.822. These values represent a sizable reduction in 
the recovery temperature Tj- and the effective temperature Tg. Thus, 

the aerodynamic heating characteristics, at high Mach numbers do not 
appear to be as severe as is indicated by simpler theories. 

The equations defining dj. and 3g are of the same form. The 

only difference between the two expressions is that the parameters 
occurring in 3g correspond to the condition of Tg* = Constant, 
whereas the values in the expression for 3^ correspond to the condi- 
tion of zero heat transfer. The functions dj. and 3g became identical 

for the condition of zero heat transfer. Since the temperature dif- 
ference Tg - Tg eirises naturally in equation (2l) for the heat-transfer 

rate, it appears that the proper reference temperature to be used for 
heat-transfer characteristics is Tg rather than T^,. 

With the heat-transfer coefficients given by equations (22) and (23), 
two Nusselt numbers may be defined as 
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Nue = ^ = -ePTi ^ 


hx r— Te* - Tg* 

NU = — = -gPTn ^ 

1 ■‘■r ■‘■s 


„ V - Tg^ 

= Niu 

T * - Tg* 


m 


Figure 6 shows that the quantity |Hu| increases with increasing Mach 
number, whereas |N\ig| decreases. The Nusselt numher |HUg) decreases 


with increasing Mach number in direct proportion to the dFunction g. 
The quantity |Nu 1 increases with Mach numher, however,' since the quan- 
Tf.* - Tg* 

tity — = ^ increases faster than g decreases. There is little 

m * _ m * ° 

■‘■r -‘■s 

change in at Mach nimibers above 7^ since the heat- transfer coef- 

ficient hg changes vei*y little. ’ The change of Nu or NUg over 
the entire Mach number range is not large. 


From equations (19) and (2l) a sing)le relation between the heat- 
transfer rate and the shear stress is found: 


q 





This relation, together with the heat-transfer rate from equation (2l), 
is shown in figure 7. Since the skin friction decreases and Jqj increases 
with Mach niimber, the ratio (qAsl increases with Mach number for a given 
value of Tg*. The calculations at a Mach number of 7 Indicate that the 
effect of increasing values of Tg* is to decrease both the heat-transfer 
rate |q| and the shear stress Tg, and the ratio lq/rg| also decreases* 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , July 19, 1951 
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Figure 2. - Velocity and temperature profiles for = 7, 
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Figure 3»- Velocity and temperature profiles for the condition of zero 

heat transfer at the surface. 
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